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ABSTRACT 

Recently, the CO J=3-2 observational result of the envelope of the 21 /im 
PPN IRAS 07134+1005 has been reported. Assuming that the CO J=3-2 line 
was optically thin, the mass-loss rate of the superwind in this PPN was found 
to be at least 2 orders of magnitude lower than the typical range. In order to 
obtain a more accurate mass-loss rate, we reexamine this data and construct 
a radiative transfer model to compare with the data. Also, in order to better 
resolve the superwind, we adopt a different weighting on the data to obtain maps 
at higher resolution. Our result shows that the CO J=3-2 emission is located 
slightly further away from the central source than the mid-IR emission, probably 
because that the material is cooler in the outer part and thus better traced by 
the CO emission. At lower resolution, however, the CO emission appeared to be 
spatially coincident with the mid-IR emission. Our model has two components, 
an inner ellipsoidal shell-like superwind with an equatorial density enhancement 
and an outer spheroidal AGB wind. The thick torus in previous model could 
be considered as the dense equatorial part of our ellipsoidal superwind. With 
radiative transfer, our model reproduces more observed features than previous 
model and obtains an averaged superwind mass-loss rate of ~ 1.8 x 10~ 5 M® 
yr _1 , which is typical for a superwind. The mass- loss rate in the equatorial plane 
is 3 x 10~ 5 M Q yr -1 , also the same as that derived before from modeling CO 
J=l-0 emission. 

Subject headings: (stars:) circumstellar matter — planetary nebulae: general - 
stars: AGB and post- AGB — stars: individual (IRAS 07134+1005) — stars: 
mass-loss 
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Introduction 



Proto-planetary nebula (PPN) is an object in a transient phase between an asymptotic 
giant branch (AGB) phase and a planetary nebula phase in the end stage of a low to inter- 
mediate mass (1— 8M Q ) star. It consists of a post-AGB stellar core (star) and an extensive 
circumstellar envelope of dust and gas. It is bright at infrared (IR) wavelength because 
of the dense and warm circumstellar dust. Recently, a large number of observations have 
revealed that PPNe are mostly bipolar, multipolar or elliptical with a torus-like structure in 
the equator. Moreover, most of them have an extensive and roughly round halo, produced 
by the progenitor AGB star (see, e.g., review by Balick & Frank 2002), indicating that there 
is a structural change from AGB to PPN phase. The shaping mechanism of PPNe is still 
unclear and believed to be closely related to the mass-loss process of the post-AGB star. 

IR AS 07134+1005 (hereafter 107134) is a well-studied carbon-rich PPN with a 21 /mi 
feature (IKwok et al.lll989l : IVolk et al.lll999l ) that is also seen in some other carbon-rich PPNe. 
The chemical study of 107134 reveals a metal-poor central star with abundant C, N, O, and 
s-process elemen ts, indicating that the star has g one through th i rd dre dge-up when it was on 
the AGB phase (IVan Winckel fc Reyniersll2000l ). iBarthes et al.l (120001 ) studied the pulsation 
of the central star revealing a variation of 36.8 days, indicating a stellar mass of 0.6 M Q . 



In IR images (IKwok et al.l 120021 ; lUeta et al.ll2005l ). 107134 was seen with an elliptical 
dusty shell with an equatorial density enhancement (i.e. a torus-like structure) surround- 
ing its central star. In the optical, the elliptical s hell was seen sur rounded by a roughly 
round halo produced by the AGB wind in the past (jUeta et al.ll2000l ). The CO J=l-0 map 
(IMeixner et al.ll2004L hereafter M2004) also showed a torus-like structure similar to that seen 
in the mid-IR image. A radiative transfer model made by M2004 indicates that the envelope 
of 107134 has two components, an inner warm and dense superwind corresponding to the 
elliptical shell, and an outer cool and sparse AGB wind corresponding to the round halo. 



Nakashima et al.l (120091 . hereafter N2009) have repo rted the CO J= 3-2 observation of 
this PPN obtained with the Submillimeter Array [SMA JHo et all 120041 )] . Unlike CO J=l- 



0, CO J=3-2 traces mainly the superwind component because it traces warmer and denser 
material than CO J=l-0. N2009 have proposed a morpho-kinematics model to compare with 
the observation. They assumed the CO J=3-2 line is optically thin and estimated a lower 
limit of the mass-loss rate, which is at least two orders of magnitude lower than the typical 
range (10~ 7 to 10~ 4 M yr _1 , see, e.g., the review by Van Winckel 2003). For studying the 
shaping mechanism of 107134, we need to know the mass-loss rate, spatial structure and 
kinematics of the superwind component as accurate as possible. Therefore, in this paper, 
we reexamine the CO data and construct a radiative transfer model to compare with the 
data. The details about our data reduction and mapping are described in Section 2. The 
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observation results are presented in Section 3. Our model is described and compared with the 
observation in Section 4. We discuss and conclude our work in Sections 5 and 6, respectively. 



2. Observation 

The SMA CO J=3-2 observation of 107134 in this paper has been reported by N2009 
and please refer to their paper for the details. In this Section, we only summarize the 
parameters of this observation in Table [U and describe the differences between our and their 
data reduction and mapping. 

8 antennas were used in the observation. We adopted the data from 7 (one more 
than N2009) of them in data reduction and mapping because antenna 7 was excluded due 
to its weak fringes and amplitudes over the whole observation. In addition, we removed 
the first half data of antenna 4 because of its scattering phase. In mapping, a "super- 
uniform" weighting is used to achieve the best compromise between the sensitivity and 
angular resolution. It results in an angular resolution of 1.66" x 1.46" with a position angle 
(PA) of 17°, which is higher by 40% in beam area than that in N2009, showing the compact 
structure clearer. Our maps will show similar features to those of N2009 if we convolve our 
maps with their angular resolution. The channel maps have a resolution of ~ 0.7 km s _1 
per channel, with a rms noise (hereafter a) of 0.27 Jy Beam -1 . They are used to produce 
integrated intensity map, position velocity (PV) diagrams, and spectrum. The channel maps 
we shown in this paper are binned to have a lower resolution of ~1.4 km s _1 per channel 
to show how the structure changes with velocity. Although this velocity resolution is lower 
by 40% than that of N2009, it is enough to show the main structure of 107134 in different 
velocity. 



3. Observational Results 



3.1. Integrated intensity map 



In Figure [Tl we s uperposed our SMA CO J=3-2 (int egrated intensity) map on the near- 
IR JUeta et al.ll2005h . CO J=l-0 (M2004), and mid-IR fcwok et alJbooj ) maps as well as 
the CO J=3-2 map of N2009 to analyze the structure of 107134. The near-IR emission 
mostly traces the light of the central star scattered by the surrounding dust, and the mid-IR 
emission mostly traces the thermal dust emis sion. The near-I R map shows a well-defined 
elliptical shell with a major axis at PA ~ 25° ( jUeta et al.ll2005f ). and the mid-IR map shows 
a roughly round envelope surrounding the central star. The maps all show a double-peak 
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structure around the minor axis (PA ~ 115°), indicating that the circumstellar envelope of 
107134 has an equatorial density enhancement (i.e., a torus-like structure). At low resolution, 
the two pea ks of the CO J=3-2 em ission appeared to be spatially coincident with those in 
the mid-IR (INakashima et al.ll2009l ). But here at higher resolution, the two peaks of the 
CO emission are located slightly further away from the source than those in the mid-IR, 
especially for the western peak. The two peaks are asymmetric about the central star, with 
the eastern peak closer to the central star than the western peak by about 20%. Note that 
both the CO J=3-2 and the CO J=l-0 emission are weak in the north of the central star 
than that in the south, in contrast to that seen in the mid-IR. 



3.2. Channel maps 

Our channel maps (Figure [2]) are similar to those of N2009 (their Figure 3), showing 
similar variation of the morphology with the velocity. For example, the size of the emission 
gradually increases from high blue- and red-shifted velocities toward the systemic velocity 
(72 km s _1 , as found in our model described later), single peaks are seen at high blue- and 
red-shifted velocities, an opening (less emission) toward the south (i.e., inverted U-shape, 
"convex upward" in N2009) from 66.4 to 69.2 km s _1 , an opening toward the north (i.e., 
U-shape, "convex downward" in N2009) from 73.4 to 79 km s -1 , and two elongated and 
clumpy peaks in the east and west near the systemic velocity (from 70.6 to 72 km s _1 ). 
Notice that, the high velocity emission peaks are not exactly at the map center (central star 
position). The highest blue-shifted peak is shifted slightly to the north, and the highest 
red-shifted peak is shifted slightly to the south-west. 



3.3. PV diagrams 

The PV diagrams cut along the major (the poles, PA= 15°, as found in our model 
described later) and the minor (the equator, PA= 105°) axes of the ellipse are presented 
in Figure [3 ] Not e that the major and minor axes here are different from those found in 



Ueta et al. 



(120051 ). These PV diagrams are similar to those of N2009 that cut along the 
similar PAs (0° and 90°, their Figure 4). At higher angular resolution, however, our PV 
diagrams show a clearer ring-like structure [as indicated by a white dashed ellipse on Figure 
[3] (b)] in both cuts. The lack of emission inside the ring-like structure is consistent with 
a detached envelope around the central star. In the PV diagrams, there are two obvious 
asymmetric emission distributions. One is that the high red-shifted emission is stronger 
than the high blue-shifted emission toward the star position in both PV diagrams. The 
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other is that in the PV diagram cut along the major axis, the south-western emission in the 
red-shifted part is much stronger than that in the blue-shifted part, producing a gap (Gap 
1) there. There is another gap near the highest blue-shifted velocity (Gap 2), separating 
the highest blue-shifted emission from the main emission structure. Note that in the PV 
diagram cut along the major axis, the emission at the highest blue-shifted and red-shifted 
emission are not at the source position, consistent with that seen in the channel maps as 
mentioned above. 

From the observation, we suggest that the envelope of 107134 is a radially expanding 
ellipsoidal shell, with the red-shifted emission stronger than the blue-shifted emission. This 
envelope is inclined to the plane of the sky, producing the U- and inverted U-shape structures 
in the channel maps. In the following section, we introduce our radiative transfer model and 
compare it with this observation in order to obtain more accurate mass- loss rate, morphology 
and kinematics of the envelope. 



4. An Expanding ellipsoidal model 
4.1. Model 

Here we construct a code with a model different from that of N2009 and include the 
radiative transfer to calculate the CO J=3-2 emission, rather than assuming an optically 
thin emission. As in M2004, the envelope has two components, an outer spherical AGB 
wind and an inner ellipsoidal shell-like superwind (Figure H]). In the superwind component, 
we do not adopt a spheroidal shell as in M2004 nor a thick torus as in N2009, because the 
near-IR image [Figured] (a)] shows an elliptical shell-like morphology in the inner part of this 
envelope. Thus, the superwind is assumed to be an ellipsoidal shell elongated in the north- 
south direction and looks like a football. In the cartesian coordinate system, the ellipsoidal 
shell can be described with the following equation, 

x 2 y 2 z 2 

with the x — y plane being the equatorial plane and the z— axis being the major axis in the 
north-south direction. Here, a and b are the one-half of the minor and major axes of the 
ellipsoidal shell with an constant ellipticity defined as e = 1 — | . The ellipsoidal shell has 
an inner radius R[ n and a thickness Aa in the minor axis, and thus a = Ri n to R m + Aa 
(which is the outer radius). On the other hand, the spherical AGB wind has an outer radius 
of .Rout; and its inner boundary is the outer boundary of the superwind. 



We assume that the envelope is expanding radially, and the expansion velocity of the 
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superwind follows the same equation as the ellipsoidal shell for self-consistency, 

K? K 2 K 2 ' l) 

with \4 and 14 being the expansion velocities in the minor (equator) and major (pole) axes, 
respectively. As can be seen from the equation, the expansion velocity increases from the 
equat or to the poles, qualita tively similar to that expected in some of the ISW simulations 
(e.g., iDwarkadas et al.lll996l ). As for the AGB wind, it expands radially with a constant 



velocity of V a . 

Below are our density and temperature profiles of the molecular gas (molecular hydro- 
gen) in the envelope (in spherical coordinate system). 



p H (R 9) = \ ^k^ 6 ^ Su P erwind component, 
1 ^ 47r#*Va AGB wind component, 



^gas (R) 



Tm Superwind component, 

Ti n F T (^^j AGB wind component, 



where M eq is the mass-loss rate of the superwind component in the equator. The density 
decreases with radius R with a power-law index of -2, in both the superwind and AGB 
wind components, and it has a sudden drop with a factor F at the superwind-AGB wind 
boundary. In order to reproduce the observed equatorial density enhancement, the density 
of the superwind component is multiplied by a simple torus function 

/(0) = l-acos0, (3) 

where 9 is the angle from the z-axis, and the value of a is between and 1. Here f(9) equals 
to 1 at the equator and 1-a at the poles. T in is the temperature at i?j n . The temperature is 
assumed to decrease with R with a power-law index of (3 in the superwind and 7 in the AGB 
wind. It also has a sudden drop with a factor Ft in the superwind-AGB wind boundary. In 
the AGB wind, the temperature is set to 5 K when it decreases to below 5 K. The density 
and temperature profiles along the minor axis of the two components are shown in Figure |5j 

Radiative transfer is used to derive the CO J=3-2 emission in our model. Thermal 
line width V% n and the line width due to the turbulence velocity V tur b are also included. 
The systemic velocity V^ ys is assumed to be a free parameter. Also, we rotate our model 
counterclockwise by a position angle (PA) from the z-axis and tilt it with an inclination angle 
i with the north part tilted away from us. To compare the model results with the observation 



properly, we use MIRIAD (ISault et al.lll995l ) to derive the model visibility from the model 
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data cubes with the observed uv— coverage, and then use the same imaging procedure as we 
did for the observed channel maps to derive the model channel maps. As a result, the model 
and the observed channel maps have the same spatial and velocity resolutions. With the 
channel maps, we can obtain the integrated map, PV diagrams, and spectrum of our model. 



4.2. Model Results 

The best-fit parameters are listed in Table [2j There are 18 parameters in our model, 
3 are constant and 15 are free. The 3 constant parameters are distance of the source, CO 
abundance (CO/H 2 in number density) and -R ou t- For the CO abundance, we adopt the 
value of 9.2 x 1(T 4 from M2004. If we adopt the value of 7.4 x 1(T 4 as in N2009, which is 
lower by ~ 20% than that in M2004, our mass-loss rate will increase only by ~ 20%. The 
15 free parameters are R{ n , Aa, e, PA, i, V sys , V a , Vt ur b, T in , Ft, (3, 7, F, M eq and a. We 
compare the results of our model with the observation, and determine the best fit parameters 
and their error bars based on the following five criteria. (1) The total flux of our model can 
not be more or less by 10% than that of the observation. Here the total flux is derived from 
within a 4" wide box centered at the central star, where the emission mainly comes from 
the superwind region. (2) The blue-shifted peak of the spectrum of our model can not be 
more or less by 10% than that of the observation. (3) The highest contour in the moment 
map of our model can not be more or less by one contour level than that of the observation. 
(4) The lowest contour in the moment map of our model should be in between the lowest 
and penultimate-low contours in the moment map of the observation. (5) The morphology 
trend in the channel maps and PV diagrams in our model should be similar to that in the 
observation as mentioned in Section 3. 

As shown in Figures |6] and [71 our model can reproduce most of the observed structures 
and kinematics of the CO envelope. For example, how the morphology and intensity of the 
emission change with velocity in the channel maps, the equatorial enhanced emission in the 
integrated intensity map, the intensity ratio of high red-shifted to high blue-shifted peaks in 
the spectrum, and the emission gaps in the blue-shifted part of the PV-diagrams. The mean 
optical depth of the CO J=3-2 emission estimated from the intensity ratio of the two peaks 
in our model spectrum toward the source position is ~1. 



In our models, the values of f?i n and Aa are very similar to those estimated by lUeta et al. 



( 120051 ) from the near-IR image, which are 1.4" and 0.8", respectively. Also, e ~ 0.2, almost 



the same as the value of 0.21 estimated from the optical image with the filter of F547M 



( jUeta et al.ll2000f ). The mean expansion velocity of the superwind is also the same as the 



expansion velocity estimated by M2004 (10.5 km s 1 ). In our model, V tnT h ~ 1 km s 1 , and 
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it is much important than K n , which has a m aximum value of ~ 0.2 km s _1 . Our PA is ~ 
15° ±5°, in between those of lUeta et al.l ( 120051 ) a nd N2009. The in clination angle is estimated 
to be ~ 25 ± 5°, in between those estimated by lUeta et al.l (120051 ) and N2009. M eq obtained 
from our model is the same as that in M2004. In our model for CO J=3-2, the superwind 
has a P of -0.8, much shallo wer than -2.0 found in CO J=l-0 (M2004) , but steeper than -0.4 
found in IR dust emission ( iCollison fc Fixlll99ll ; iMeixner et al.l 119971 ). On the other hand, 
in our model for CO J=3-2, the AGB wind has a 7 of -1.5, much steeper than -0.25 found 
in CO J=l-0 (M2004). We will discuss the differences later in the Discussion Section. 

We can then derive 5 other quantities from our model: dynamical time at the inner 
boundary of the superwind (tdyn)j duration of the superwind (£ sw ), averaged mass-loss rate 
of the superwind (M sw ) averaged over 9, as well as duration (t ag b) and mass-loss rate (Magb) 
of the AGB wind. We divide R in by V a to obtain t dyn of ~1590 years, divide Aa by V a to 
obtain t sw of ~1100 years. Our Msw is 1.8 x 10~ 5 M© yr _1 , which is lower by 40% than 
that of M2004. We divide the thickness of the AGB wind (~5.5") by V a to obtain £ agb of 
~ 6760 years. Our tdyn is slightly larger than that obtained by M2004 (1240 years), with a 
larger inner radius than theirs. Our t sw and t ag b are similar to those estimated by M2004, 
which are 840 and 6570 years, respectively. We multiply M eq by F to obtain Magb of 3.6 
(± 1.5) x 10~ 6 M yr _1 . Magb derived by M2004 is 5.1 x 10~ 6 M yr" 1 and it equals our 
upper limit. In our model, the CO emission of the AGB wind is pretty weak, with less than 
3a even without having applied the observed uv coverage, and thus will not be detected in 
our observation. 



On the other hand, our model can not reproduce the following observed features: 



1. Our model is symmetric and thus can not reproduce the different radii between the 
eastern and western peaks in the observed integrated intensity map. The reason for 
this asymmetric is unclear, and will be discussed later. 

2. The intensity excess at low red-shifted velocity in the observed spectrum [Figure [3(b)]. 
The spectrum is obtained by averaging over a circular region of 1" in diameter around 
the center, showing mainly the emission from the front and back walls of the envelope. 
Thus, the intensity excess indicates that there could be an additional material with 
low expansion velocity in the back wall (red-shifted part) of this envelope. 

3. In our model, the highest blue-shifted emission is in the southwest and highest red- 
shifted emission is in the northeast of the center (Figure E]), in opposite to those seen 
in the observation. The superwind could have a non-radial velocity component, for 
example, a poloidal velocity (directed from the equator to the poles). In the front wall 
of the superwind envelope, a poloidal velocity could increase the projected velocity of 
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the north-eastern gas and decrease the projected velocity of the south-western gas near 
the map center. Same thing is for the back wall of the superwind envelope but with 
an opposite velocity sense. Alternatively, there could be an additional high-velocity 
component not included in our model. 



5. Discussion 

5.1. Superwind material in IR and CO observations 

As compared to N2009, the superwind is better resolved in our map optimized at a 
higher angular resolution. As a result, we can now study the superwind in more details. At 
a higher rotational transition, CO J=3-2 can trace warmer and denser material than CO 
J=l-0. The refore, our CO J=3-2 map can reveal the superwind material better than the CO 



J=l-0 map (IMeixner et al.ll2004j ). which shows both the superwind and AGB wind materials. 
Moreover, as mentioned before in Section 3, the two emission peaks in the equator in our 
CO J=3-2 map are located slightly further away from the source than those in the mid-IR 
map. It is likely because that the mid-IR emission mainly traces the warmest material in 
the innermost part of the superwind, CO J=3-2 emission mainly traces the warm material 
in middle part of the superwind, and CO J=l-0 emission traces the cool material further 
out in the superwind and the AGB wind as well. Thus, in the superwind, the different 
temperature power-law index (3 in different tracers and transitions may suggest that the 
power-law index becomes steeper with the distance from the source. On the other hand, the 
temperature power-law index of our AGB wind (7 = —1.5) in CO J=3-2 is in betw een those 



of the superwind (-2.0) and AGB wind (-0.25) in CO J=l-0 (IMeixner et al.1 120041 ). This is 
probably because that the superwind and AGB wind in the CO J=l-0 observation can not 
be clearly separated at low resolution. It is also possible that our observation is not sensitive 
to the AGB wind. 



5.2. Comparison with N2009 results 



N2009 has proposed a model to compare with this observation. Their model contains 
an inner toroidal superwind and an outer spheroidal AGB wind. Figure H] right shows their 
model on top of our model. It is clear from the figure that their torus could actually be 
considered as the dense equatorial part of our ellipsoidal superwind. Note that the radial 
extent of their torus is actually about one time larger than our superwi nd, and the PA der ived 
by our model is between that derived by N2009 and that estimated by lUeta et al.l (120051 ). In 
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their model, there was no radiative transfer and they convolved their channel maps with a 
circular beam rather than the observed elliptical beam. Their model reproduced the U- and 
inverted U-shaped morphologies in the channel maps and the ring-like PV structure in the 
PV-diagram. However, in order to reproduce the two elongated and clumpy structures in 
the equator around the systemic velocity, their model needed a distorted torus. Moreover, 
unlike that seen in the observation, the emission intensity in their channel maps did not 
gradually increase from the blue-shifted to red-shifted channels, but had a sudden increase 
in some blue-shifted channels. In their PV-diagrams, Gap 2 was not seen on the high blue- 
shifted velocity side. Our model, on the other hand, can reproduce all the above observed 
features reasonably well. Furthermore, M eq estimated from our model is the same as that 
in M2004 and is in the typical range for an AGB star. In contrast, N2009, by assuming 
that CO J=3-2 line is optically thin, estimated a lower-limit mass- loss rate of ~ 10~ 9 M 
yr _1 , which is 4 orders of magnitude lower than our M eq derived from our radiative transfer 
model. Therefore, radiative transfer is really needed to estimate the mass-loss rate properly. 
Note that, however, it is unclear how N2009 obtained such a low mass-loss rate since the 
optical depth toward the central position in our model is only ~ 1. 



5.3. Shaping mechanism of elliptical PPN 



As mentioned in the Introduction Section, most of the low- to intermediate-mass stars 
experience a structural change fro m spherical to ell i ptical or to bipola r structure when they 
evolve from AGB to PPN phase (jUeta et al.l |2000| ; ISahai et al.l 120071 ). It is still uncertain 
what mechanism can cause this structural change. Moreover, most of PPNe contain an 
additional density enhancement (torus-like) structure that is perpendicular to their elliptical 
or bipolar structure. 107134 is one of the typical elliptical PPNe. It has an elliptical shell 
with an equatorial density enhancement embedded in a round halo, indicating that it also 
experiences the structural change mentioned above. At the present day, one of the popular 
models to explain the forming mechanism of the density enhancement is a model with a 
binary system. In such a model, if the separation of the binary, wind velocity ejecting from 
the primary star and the mass of the companion are in appropriate range, the material 
ejecting from the primary star can accumulate around the equator with higher expansion 
velocity than the p oles, forming an asymmetric d ensity enhancement in the equator as seen 
in the observation (IMastrodemos fc Morris! 119991 ). However, in this scenario, the elliptical 
shell will appear oblate with the long axes lying in the orbital plane, in contrast to that seen 
in 107134. 



The elliptical structure of 107134 suggests that the poles of the superwind have a higher 
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expansion velocity than the equator. Thus if the binary system is to work, an additional 
element will be needed to increase th e expansion velocity at the po l es. A bipolar jet launched 
by accretion effect of the companion (IMastrodemos fc Morrislll998l ; iFrank &: Blackmanll2004l ) 
is a candidate that can produce high expansion velocity at the poles and it h as been used to 
explain the shaping mechanism of bipolar PPNe (see, e.g.. lLee fc Sahaill2003l ). However, the 
observations toward 107134 do not show any high-velocity emission. It is possible that the 
jet could be an atomic jet and thus could not be detected in CO observation. Alternatively, 
a tenuous (and thus unseen) isotropic post-AGB wind with higher expansion velocity than 
the superwind coul d have been launched a nd interacted with the superwind, producing the 
elliptical structure (IDwarkadas et al.lll996l ) of 107134. In either cases, 107134 could be in the 
transient phase changing from a round to a bipolar structure. 



However, the long-term (20 years) radial velocity observation of iHrivnak et al.l (1201 if ) 
seems to not support the binary scenario in 107134. If this is the cas e, the structural c hange 
of 107134 could be caused by some mechanisms due to a single star. iDorfi fc Hoefnerl (119961 ) 
proposed that a low rotation of a single AGB star would induce a preferential mass loss 
with higher velocities in the equatorial plane. However, in this scenario, th e envel ope would 
form an oblate shape, inconsistent with that seen in 107134. iMatt et al.l ( 120001 ) proposed 
the dipole m agnetic field in an AGB star could lead the mass loss along the equator. But 
Sokerl (120061 ) argued that the magnetic field would carry angular momentum away from the 
stellar envelope faster than the mass is lost by the wind, causing the star to spin down on 
a short timescale. In this case, the lifetime of magnetic field might not be long enough to 
produce the equatorial density enhancement structure. ISokerl (119981 . 120001 ) proposed that 
the concentration of cool magnetic spots toward the equator on the surface of an AGB star 
would lead to the equatorial density enhancement, because dust would form above the cool 
spots more than other area. Moreover, the mass loss velocity above the cool spots is lower 
than other directions due to the weak radiation from the cool spots, which has potential to 
form an elliptical envelope like 107134. However, it is unclear how above models can explain 
the asymmetric density distribution in the envelope of 107134. 



6. Conclusion 

We have reexamined and remodeled the SMA CO J=3-2 observational result of 107134. 
Our main results are the following. 

1. The CO map is consistent with the IR maps, showing an equatorial density enhance- 
ment. The CO emission is located slightly further away from the central source than 
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the mid-IR emission, probably because that the material is cooler in the outer part 
and thus better traced by the CO emission. 

2. Our model has two components, an inner ellipsoidal shell-like superwind with an equa- 
torial density enhancement and an outer spheroidal AGB wind. The thick torus pro- 
posed by N2009 could actually be considered as the dense equatorial part of our ellip- 
soidal superwind. 

3. Our model can reproduce the observation reasonably well, better than the model pro- 
posed in N2009. The superwind mass-loss rate in the equator is estimated to be ~ 3 
x 10" 5 M Q yr" 1 , the same as that derived by M2004 and is in the typical range for an 
AGB star. The mean expansion velocity of the superwind and AGB wind are 10.5 and 
9.3 km s _1 , respectively. The mass-loss durations of the superwind and AGB wind are 
1100 and 6760 years, respectively, similar to those estimated in M2004. The superwind 
ended its ejection about 1590 years ago. 

One of the popular models to explain the forming mechanism of the density enhancement 
is a model with a binary system. However, in order for this binary model to work for 
this PPN, an additional element, such as a bipolar jet or a tenuous post-AGB wind, may 
be needed to produce the prolate elliptical structure with the density enhancement in the 
equator. Alternatively, the forming mechanism of the density enhancement could be due to 
a single star. 
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Fig. 1. — Integrated intensity map of the S MA CO J=3-2 em ission (green contours) super- 
posed on (left) the n ear-IR (1.1/xm) image ( lUeta et al.l 120051 ). (center) the CO J=l-0 map 
( iMeixner et al.ll2004j ). (right) the integrated CO J=3-2 int ensity map of N20 09 (white con- 
tours) and the mid-IR (12.5/mi) image [color background (IKwok et al.l 120021 )]. respectively, 
in the same scale. The integrated intensity map is integrated from 57.6 to 86.5 km s _1 . The 
"+" symbol marks the position of the central star. The first and last contours of the CO 
J=3-2 emission are 2.8a and 66.2<r, respectively, and the middle contours are from 7a to 63cr 
with a step of 7a, where la is 1.23 Jy beam -1 km s _1 . The first contour level is the same 
as that of the integrated map of N2009. The resolutions of the CO J=3-2 emission of this 
work (1.66" x 1.46", PA= 17°) and of N2009 are shown in the bottom-right corner of (left) 
and (right), respectively, and the resolution of CO J=l-0 emission is shown in bottom- left 
corner of (center). 
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Table 1. SMA observations of IRAS 07134+1005 



Parameter 

R.A. (J2000.0) a 
Decl. (J2000.0) a 
Primary beam (FWHM) 
Velocity resolution (AV) 
Flux calibrator 
Passband calibrator 
gain calibrators 
Synthesized beam 
a (AV ~ 0.7 km s" 1 ) 



Number 



07 h 16 m 10?26 
+9°59'48"0 
~ 34" 
~ 0.7 km s _1 
Callisto 
Callisto 
0739+016, 0750+125 
1.66" x 1.46" (PA = 17°) 
0.27 Jy Beam" 1 



a T he central star is assu med to be at this coordinate, 
as in iMeixner et al.l ( 12004] ) 
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Table 2. Model Parameters 



Parameter 


Value 


Reference 


Constant 


Distance 


2.4 kpc 


1 


CO/H 2 


9.2 xKT 4 


2 


-Rout 


~ 19300 AU (8") 


2 


Free parameters 


-Rin 


1.3" ± 0.1" (~3100 AU) 




Aa 


0.9" ± 0.1" (~2160 AU) 






0.2 ± 0.05 




PA 


15° ± 5° 




i 


25° ± 5° 




^sys 


72 ± 0.35 km s" 1 






9.3 ± 0.35 km s" 1 




Kurb 


1 km s _1 




T- 

J - in 


70 ± 5 K 




pip 


0.25 ± 0.1 




P 


-0.8 ± 0.2 




7 


-1.5 ± 0.5 




F 


0.12 ± 0.05 




Meq 


3 (± 1) xlO" 5 M yr- 1 




a 


0.8 ± 0.1 




Output 


^dyn 


~ 1590 yr 




tsw 


~ 1100 yr 




^agb 


~ 6760 yr 




averaged Msw 


1.8 (± 0.6) x 10~ 5 M yr" 1 




Magb 


3.6 (± 1.5) x lO" 6 M yr" 1 





Ref erences. — (1) iKnapp et all ( 120001 ); (2) Meixner et al 
j2004h . 
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Fig. 2. — Observed channel maps in contours and grey scale. The contour levels are from 
4.3cr to 39a with a step of 4.3a, where a is 0.19 Jy Beam -1 . Top-left corner of each channel 
map shows the LSR velocity in km s -1 . The systemic velocity is in the channel of 72 km s _1 . 
The velocity width is ~ 1.4 km s -1 . The white cross symbol at the center of each channel 
map marks the position of the central star. The resolution is shown in the bottom-left corner 
of the first channel map. The wedge shows the intensity scale in Jy Beam -1 . 
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Fig. 3. — Observed PV diagrams cut along (a) PA= 15° and (b) PA = 105°, respectively, 
in contours and grey scale. The contour levels are from 3a to 23<r with a step of 3a, where 
a is 0.27 Jy Beam -1 . The capital NE, SW, SE and NW indicate the north-east, south-west, 
south-east and north-west directions on the observed integrated intensity map. The wedge 
shows the intensity scale in Jy Beam -1 . 
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Fig. 4. — Left: A transect map of the density (grey scale) and velocity (grey arrows) distri- 
butions of the molecular envelope in our model (Section 4). This map is taken from a plane 
along the symmetric axis (PA=15°). The inner dense elliptical region is the superwind com- 
ponent and the outer sparse spherical region is the AGB wind component. Right: A transect 
map of the model of N2009 on top of ours in the same scale but taken from a different PA, 
theirs is 8° and ours is 15°. The grey circle shows the outer boundary of their spheroidal 
AGB wind and the two hatched rectangular regions show their thick torus-like superwind. 
In the right sketch, the axes and labels are for the model of N2009. 
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Fig. 5. — Equatorial density (upper panel) and temperature (lower panel) of the molecular 
function of radius for the superwind and AGB wind components in our model. 
Please note that we do not show the entire AGB wind region, which has an outer radius of 
~ 3 x 10 17 cm. 
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Fig. 6. — The channel maps derived from our model. The contour levels, resolution, velocity 
width and other map information are the same as those in Figure [2j 
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Fig. 7. — Comparison of our model and the observation in the integrated intensity map, 
spectrum, and PV diagrams. Grey contours and spectrum are from the observation. Orange 
contours and spectrum are from our model, (a) shows the integrated intensity maps. The 
contour levels are the same as those in Figure [TJ Note that however, we do not include here 
the last contour of Figure HJ (b) shows the spectra toward the central star position averaged 
over a circular region with a diameter of 1". (c) and (d) show the PV diagrams cut along 
PA = 15° and PA = 105°, respectively. The contour levels are the same as those in Figure 

El 



